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Evidence is presented for the necessity of including duality violations in a consistent description of spectral 
function moments employed in the precision determination of a a from r decay. A physically motivated ansatz 
for duality violations in the spectral functions enables us to perform fits to spectral moments employing both 
pinched and unpinched weights. We describe our analysis strategy and provide some preliminary findings. Final 
numerical results await completion of an ongoing re-determination of the ALEPH covariance matrices incorpo- 
rating correlations due to the unfolding procedure which are absent from the currently posted versions. To what 
extent this issue affects existing analyses and our own work will require further study. 
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1. Introduction 

Hadronic decays of the r lepton provide a par- 
ticularly clean environment for the study of low- 
energy QCD. The mass of the r (which is the only 
lepton heavy enough to decay into hadrons) is 
large enough that pcrturbativc QCD should pro- 
vide a good description of the inclusive decay pro- 
cess, but not so large that nonpcrturbativc effects 
can be completely neglected. 

Since the seminal work of Ref . p] , it has been 
recognised that the perturbative series supple- 
mented with the Operator Product Expansion 
(OPE) provides a solid framework for describ- 
ing these decays. The main achievement of this 
program is a precision determination of the QCD 
coupling a s . The determination of a s at differ- 
ent scales provides a highly nontrivial test of the 
theory Hadronic r decays also yield values for 
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vacuum condensates [5] , low-energy constants [3] , 
and the CKM matrix element \V US \ [4]. 

Unfolded high-statistics spectral functions ex- 
tracted from the final LEP data are available from 
ALEPH [2] and OPAL [5]. On the theory side, the 
0{a A s ) term of the perturbative contribution has 
been calculated recently [6]. This has triggered 
several reanalyses of hadronic r decays from the 
ALEPH and OPAL spectral functions. While the 
average value of a s from t data is competitive 
with the current world average, the values ob- 
tained by different groups are barely compatible. 
After the evolution to the Z mass, a sample of 
results obtained by various groups is 

0.1202 (6)exp(18)th i , 
0.1212 (5) oxp (9) th [7] , 
0.1180 (4) exp (7) th , 
0.1187 (6) ex p(15)th M ■ 
Furthermore, there is some tension between the 
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result of Ref. [7] and recent lattice determina- 
tions; e.g., the analysis of Ref. [TU] gives 



a s (m%) = 0.1184 ±0.0006. 



(2) 



Even though much of the variation in Eq. ([T]) re- 
sults from differences in the prescriptions used 
for resumming the perturbative series (CIPT ver- 
sus FOPT |HI8ll2j ). at the current level of pre- 
cision nonperturbative effects thus far neglected 
may become significant. 

One can think of several different, but not com- 
pletely independent, sources of systematic uncer- 
tainties in the theoretical description of hadronic 
r decays. Examples are the truncation in powers 
of a s and/or variation in the choice of resumma- 
tion prescription for the perturbative series |8ll2j , 
unchecked assumptions regarding the truncation 
in dimension of the OPE [5J, and contributions 
from duality violations (DVs) |13ll4j . In this 
work we address the latter. The OPE is believed 
to be an asymptotic expansion at best, and it 
breaks down on the Minkowski axis. This lack 
of convergence is related to the presence of DVs, 
and indeed, our preliminary results confirm previ- 
ous observations on the presence of nonncgligiblc 
DVs for certain weighted spectral integrals [16] . 
A full reanalysis including all systematic effects 
noted above will be presented in a forthcoming 
article. 

2. Theoretical framework 

We start from the total nonstrange branching 
ratio i?f =0 defined as 

R*=° = r[T -? hadronsS= ° V A = R v + R a , (3) 
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With vector (V) and axial-vector (A) currents 
J%{x) = u-f^x) and J A (x) = uj^j 5 d(x), the 
J = 0, 1 parts of the V and A current two-point 
functions, Hy A , are defined by 



n^(?) = i J d 4 xe^(0\T{J^ A (x)J^ A (0y}\0) 

= {<f<f - q 2 gn^ { v! A (q 2 ) + <ffn$ A tf) ■ (4) 

The corresponding spectral functions Py\ = 
— Im A can be extracted from the differential 



distribution dRY' A /dq 2 , which is experimentally 



available 
dR v < A 



ds 
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Except for the pion-pole contribution to p A 1 \s), 
the J = part of the spectral function is nu- 
merically negligible. The combination Uv,a(s) = 
ITy^ '' (s) is free of kinematic singularities and an- 
alytic in the complex s-plane cut along the posi- 
tive real axis. Cauchy's theorem thus yields the 
finite-energy sum rule (FESRjl 



so 

R^(S ) = 127T 3 J 



ds 
so 



w{s/s ) pv,a{s) 
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— w(s/s )TIv,a(s) , (6) 
so 



|«|=»0 

valid for an arbitrary analytic weight w(z). 

For so 3> Aq CD , one can expect Tly,A to be 
well approximated by the OPE, i.e., 



IL VtA (s)=Il v ^(s) + AvAs), 



(7) 



with Av,a{s) a small (but in general nonzero) cor- 
rection accounting for the presence of DVs. 

Assuming DVs to vanish sufficiently fast for 
\s\ — > 00 in the whole complex plane, one can 
show that the correction from Av,a(s) to the con- 
tour integral in Eq. © can also be written as [TS] 



12tt 



ds 
so 



w(s/s ) Im A v,a(s) . (8) 



FESRs thus provide constraints on ImA^(s) 
beyond those from the region sq < s < m~ 
obtained by fitting to the experimental spectral 
function. In standard r-decay-based determina- 
tions of a s , Av,a{s) is typically neglected, with 
at most a check on the self-consistency of this 



2 In the remainder the factor of Sew I Kid 1 2 is absorbed 
into Tly^Ais). 
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assumption [5] 3 It is our aim to quantitatively 
determine the impact of DVs on the values for the 
QCD parameters accessible through such fits. 

Since at present no theory for DVs exists, we 
adopt the following ansatz for the asymptotic 
form of the DVs as a working hypothesis |15j 



-ImAy jj4 (s) — > K V . A e lv - AS $ai(av,A+Pv,As) , 

7T large s 



(9) 



reflecting the presumed asymptotic character of 
the OPE expansion. Model studies suggest the 
exponential decay, originating from the finite 
width of the resonances. The oscillatory be- 
haviour arises naturally in a spectral function 
with resonances distributed with some periodic- 
ity as, for instance, on the daughter trajectories 
in Regge theory |13ll4ll5j . With the (eight) pa- 
rameters kv,a, lv,A, ocv,a, and Pv,A we expect 
the ansatz of Eq. (JSJ) to capture the generic fea- 
tures of the DVs, thus avoiding any more specific 
models. This ansatz is certainly more general, 
and hence more likely to be realistic, than that 
used in analyses which entirely neglect DV con- 
tributions (which correspond to the special case 

Ky = hCA = 0). 

With our description of ImA^fs), the ex- 
perimental V and A spectral functions constrain 
the DV parameters, provided the asymptotic be- 
haviour of Eq. ^) has set in for s > s TO j„ with 
Smin < m r- Above s m i n one then has 



PV,A{s) =9(8- S min ) 



+ Ky A& 



Nc 

12^2 



[i + pOO] 



1 sm(a v ,A + Pv,as) 



.(10) 



The function p(s) contains the perturbative cor- 
rections up to and including 0(oq) as well 
as condensate contributions. The latter were 
shown to be numerically irrelevant for s m j n > 
1.1 GeV 2 HI]. In Ref. [15L fits to ALEPH V 



3 An exception is Ref. [7], which discusses the relevance 
of this term, though only for the combined V + A corre- 
lator, and concludes that its contribution lies within the 
error bars. On the basis of the analysis of Ref. 1151 and 
the present study we believe this conclusion needs to be 
reconsidered. 



and A spectral functions employing Eq. (|10[) have 
been carried out, testing the onset of the asymp- 
totic behaviour of Eq. © in real data. The re- 
sults show the model description of ImAv^fs) 
to be reasonable, and certainly compatible with 
ALEPH data. Similar tests confirm its compati- 
bility with the OPAL data. 

3. Strategy of our analysis 

With the OPE plus DV representation for the 
relevant correlators, Eq. ([8]) takes the form 



ds 



R^ A (s ) = 6mf -w(s/s )Uy^(s) 

\s\=s 

+ vM(so). 



(11) 



The LH side is to be evaluated using p\f A ex- 
tracted from the LEP experimental data, while 
the RH side contains the various OPE parame- 
ters as well as the parameters of the DV ansatz. 
Because of the high quality of the spectral data, 
use of a range of sq with sq > s m i n and a set of 
weight functions Wi(x) allow fits to obtain both 
OPE and DV parameters. 

Different weight functions w(x) emphasise dif- 
ferent terms of the OPE of Lly^s) ■ For want of 
a good description of DVs, existing analyses have 
been restricted to the case of so-called pinched 
weights (polynomials in i = s/sq having zeros 
at s = so which suppress contributions to the 
RHS of Eq. ^ from the vicinity of the timelike 
point on the contour \s\ = so, where DVs are 
expected to be largest). Polynomials of higher 
degree, which provide more "pinching," are sen- 
sitive to higher order terms in the OPE. In some 
analyses, for practical reasons, the OPE has been 
truncated at dimensions below that for which 
non-a s -suppressed contributions are in principle 
present. This can introduce an uncontrolled sys- 
tematic uncertainty |9|15j . 

The three main features of our strategy are 
as follows. First, because we include DV con- 
tributions, our sum rules need not be restricted 
to pinched weights only. Low-degree unpinched 
weights, sensitive to only a small number of 
condensates (modulo as-suppressed logarithmic 
corrections) are employed together with pinched 
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ones. Second, Eqs. (|5J) and (jTU)) allow simulta- 
neous fits to both the spectral function and the 
moments i?^| 1 (so). Third, following Ref. [3], we 
do not restrict ourselves only to Sq = m^. The 
use of a range of so values allows for a better use 
of the available data, facilitating the separation 
of condensate contributions of different dimen- 
sion, and providing a consistency check on the 
fits from which a s (m^.) and the condensates are 
extracted. Values for all fitted parameters should 
be independent, within errors, of the precise set 
of So an d weight functions used in their determi- 
nation. We emphasise that this check has never 
been performed in full beforeH 

In our framework, in order to determine the 
various DV parameters, one must perform a sep- 
arate analysis of V and A channels. With the 
eight parameters of Eq. ([9]) under control, one can 
then also analyse V + A or V — A. This allows 
us to consider the first Weinberg sum rule. The 
fact that our fits satisfy this constraint provides 
a consistency check on the V and A fits. 

A comment is in order regarding an important 
technical issue. In real life, data for the spectral 
functions are available in the form of binned his- 
tograms and the spectral integrals in Eq. (J6j> are 
approximated by a finite sum over bins. Even 
in the absence of correlations in the data, two 
different moments Ry" A (si) and Ry A (s2) with 
S2 > si will be correlated since they share con- 
tributions from the bins k with s^ < si. This 
correlation is stronger if S2 and si arc closer to 
each other. In practice, if one chooses adjacent 
bins, the correlation can be well over 90%. Such 
strong correlations may generate very small eigen- 
values in the correlation matrix for the moments 
Ry A (si). This renders the task of constructing 
reliable fits and a trustworthy treatment of uncer- 
tainties rather involved. Nevertheless, strategies 



4 In Ref. [9], a window of so values was used. The restric- 
tion to so > 2 GeV 2 , however, imposed by requiring the 
neglect of DVs to be self-consistent, made it impossible 
to fit the gluon condensate (a s G 2 ), which was therefore 
taken as external input. The strong anticorrelation found 
between the fitted a B and input (a a G 2 ) creates a system- 
atic uncertainty which can only be avoided in a framework 
like ours, which includes DVs, and allows also (a a G 2 ) to 
be obtained through a fit to data. 



to deal with such difficulties exist |17|18| , and can 
be used to study the quality of the fits. 

4. Evidence for duality violations 

To illustrate our argument, let us concentrate 
on fits using the weights wi(x) = 1, w<i(x) = 1—x, 
and 11)3(2;) = (1 — x) 2 . Since w\ is unpinched, 
its moments Ry A are expected to be more sen- 
sitive to DVs. The weights W2 and W3 are singly 
and doubly pinched, respectively The results dis- 
cussed in this section come from fits using 52 
values of so ranging from 1.5125 GcV 2 up to 
2.7875 GeV 2 . In the perturbative contribution 
we employ the Fixed Order prescription for the 
RG rcsummation. 

In the fits using uii(x) the parameters are 
a s (m 2 ) and the eight DV parameters Ky,A, Jv,A, 
ay, a and fiv,A- (This moment is sensitive to 
the gluon condensate only through tiny logarith- 
mic corrections.) A comparison of the fit results 
for wi(x) with the experimental moments from 
ALEPH data is shown in Fig. [TJ The main fea- 
ture of this figure is that a model in which DVs 
are neglected (kv,a = 0, red line in Fig. [1]) can- 
not account for the data. The results obtained 
with DVs (blue line), on the other hand, com- 
pare very well with the data. For the A channel 
this is true right up to the kinematic endpoint. 
For the V channel, the agreement is less good 
above 2.8 GcV 2 . Here one should bear in mind 
that (1) high-statistics measurements of 47r modes 
from BaBar suggest the ALEPH V spectral func- 
tion may be overestimated towards the end of the 
spectrum [19], and (2) the OPE plus DV fit to the 
OPAL V channel data displays no such problem. 
This point will have to be clarified in the future. 

Results for V channel fits using the W2 and W3 
moments are shown in Fig. [2j In these fits, the 
value of the gluon condensate is fitted as well. 
Additionally, W3 allows for a fit of the D = 6 
condensate. The fits with DVs (blue lines) are in 
all cases excellent. For W2, the fit including DVs 
is superior to that without. While for w$ both 
fits are of excellent quality, we recall that the fit 
without DVs requires (a s G 2 ), which can only be 
determined in a fit which includes DVs as input. 

From such fits, values for a s , the condensates 
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Figure 1. Fits to the wi(x) = 1 spectral inte- 
grals for the V channel (top) and A channel (bot- 
tom). The blue (solid) lines show fits including 
DVs whereas the red (dashed) lines represent the 
model without DVs (ki/.a = 0). 
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Figure 2. Fits to the V channel W2 (top) and 
W3 (bottom) spectral moments. The blue (solid) 
lines show fits including DVs whereas the red 
(dashed) lines represent the model without DVs 
(/c v = 0). 



and the DV parameters can in principle be ex- 
tracted. For the time being, we have not carried 
through the final version of this analysis because 
of an inconsistency we uncovered in the correla- 
tion matrix publicly available from ALEPH [20 . 
This problem is explained in the next section. 

5. Correlations in the ALEPH spectral 
function data 

In order to study the uncertainties associated 
with our fits, a Monte Carlo generator of toy data 
sets was built based on ALEPH's spectra and cor- 
relations [20] The unfolded ALEPH data are 
quite correlated, as can be seen from the upper 
panel of Fig. [3j whose data is much less scat- 
tered than the plotted errors would naively sug- 
gest. The toy data samples exhibit points that 
are much more scattered than the original data, 

5 We thank A. Hocker for this suggestion. 



as can be seen in Fig. [31 which shows the real data 
together with a representative toy sample (both 
plotted with the same errors) . The strong correla- 
tions in the upper panel are not present in the toy 
sample, which points to a problem in the posted 
correlation matrix. According to the authors of 
the original analysis, the publicly available data 
do not contain correlations due to unfolding [2Tj . 

The effect of the missing correlations seems 
to be small for the doubly or triply pinched 
weights used in recent analyses based on the pub- 
licly available ALEPH spectral data@ More pro- 
nounced effects, however, may be expected for 
FESRs based on unpinched weights, which are re- 
quired for a reliable exploration of DVs. A precise 
quantification awaits the completion of a reanal- 
ysis of the ALEPH covariance matrices. 

6 An exploratory reanalysis by the authors of Ref. [7] based 
on one exclusive channel (t~ — > tt~tt is t ) suggests that 
this is indeed the case 1211 . 
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Figure 3. ALEPH data [20] for the vector spectral 
function (top) and a MC toy data sample based 
on the ALEPH correlation matrix (bottom) . 
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